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By numerically simulating an implementation of the quantum baker's map 



l/^ ■ on an NMR quantum computer based on the molecule trichloroethylene, we 

, demonstrate the feasibility of quantum chaos experiments on present-day 

Q I quantum computers. We give detailed descriptions of proposed experiments 

^ I that investigate (a) the rate of entropy increase due to decoherence and (b) 

■ the phenomenon of hypersensitivity to perturbation. 



I. INTRODUCTION 



The quantum baker's map [|l],0 is a simple map invented for the theoretical investiga- 
tion of quantum chaos. Its mathematical properties have been studied extensively, but no 
experimental quantum systems are known which embody it. Recently, one of us has shown 
that the quantum baker's map has a simple realization on a quantum computer In this 
^ ', paper, we present realistic numerical simulations of an NMR quantum computer using 
three quantum bits (qubits) to explore the chaotic properties of the quantum baker's map. 

In Section ^ we review the definition of the quantum baker's map and its qubit realiza- 
tion. We present a simple argument [0] showing that the quantum baker's map is equivalent 
to a shift map [^,^ on a string of qubits. This leads to the definition of a simplified quantum 
baker's map which we use in the later parts of this paper. 

In Section |TTl| we review the 3-qubit NMR quantum computer based on the molecule 



trichloroethylene used in Ref. [jTO[ and give the radio-frequency (RF) pulse sequence for 



the proton and carbon spins of the molecule implementing the quantum baker's map. Since 
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decoherence cannot be neglected for this experiment, we model the NMR system by a master 
equation of Lindblad form including the Hamiltonian time evolution, the RF pulses, 
and phase noise due to the environment, using the actual experimental parameters. 

Finally, in Section |I^ we propose two specific quantum chaos experiments. In both 
experiments, we compare the quantum baker's map with a trivial regular map. One experi- 



ment analyses the rate of increase of the von Neumann entropy due to decoherence |12,13 



which is related to a quantum generalization of the Kolmogorov-Sinai (KS) entropy P,P|. 
The other experiment examines whether the 3-qubit quantum baker's map is hypersensitive 
to perturbation. Hypersensitivity to perturbation is an information-theoretic criterion for 
classical and quantum chaos [p!5|-[T8| which has been shown to be equivalent to a standard 
definition of classical chaos under general assumptions . 



II. THE QUANTUM BAKER'S MAP AS A SHIFT MAP ON QUBITS 



The classical baker's transformation [jl9|, which maps the unit square < g,p < 1 onto 



itself, has a simple description in terms of its symbolic dynamics ^0\. Each point in phase 



space is represented by a symbolic string s = ■ ■ ■ s_2S_iSo.SiS2 ■ ■ ■ where = or 1. In 
the string s, the bits to the right of the dot are the binary expansion of the q coordinate, 
and the bits to the left of the dot, read backwards, are the binary expansion of the p 
coordinate. Written formally, s is identified with a point {q,p) in the unit square by setting 
q = Sfc2~^ and p = J2k'=o^-k'^~^~^- The action of the baker's map on a symbolic string 
is given by the shift map U defined by {Us)k = Sk+i, which means that, at each time step, 
the entire string is shifted one place to the left while the dot remains fixed. Geometrically, 
if q labels the horizontal direction and p labels the vertical, the baker's map on the unit 
square is equivalent to stretching the q direction and squeezing the p direction each by a 
factor of two, then stacking the right half on top of the left. The definition of the baker's 
map through its symbolic dynamics emphasizes its prototypical character for investigations 
of chaotic maps: A very large class of chaotic maps can be shown to be equivalent to shifts 



on symbolic strings [^. Furthermore, it is shown below that the quantum baker's map is 
equivalent to a shift on a string of qubits. 

To define the quantum baker's map we quantize the unit square as in P,pT|. To 
represent the unit square in D- dimensional Hilbert space, we start with unitary "displace- 
ment" operators If and V, which produce displacements in the "momentum" and "position" 



directions, respectively, and which obey the commutation relation 



UV = VUe , (1) 

where = 1. We choose e = ^"^"^1^ . We further assume that = = 1, i.e., periodic 
boundary conditions. It follows [0,0] that the operators U and V can be written as 

f/ = e2"*^ and V = e-^^^P _ (2) 

The "position" and "momentum" operators q and p both have eigenvalues j/D, j = 
0,. 

In the following, we restrict the discussion to the case D = 2^, i.e., the dimension of 
Hilbert space is a power of two. For consistency of units, we let the quantum scale on "phase 
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space" be 27rh = 1/D = 2 ^. A transformation between the position basis {ji^j)} and the 
momentum basis is effected by the discrete Fourier transform Fjy, defined by 



fc=0 V k=0 

The D = 2^ dimensional Hilbert space modehng the unit square can be reahzed as the 
product space of qubits (i.e. N two-state systems) in such a way that 

\qj) = |aiv_i) (g) \aN-2) ® ■ ■ • ® |ao) , (4) 

where j = J2 dk'^^, cik G {0, 1} (A; = 0, . . . A^— 1), and where each qubit has basis states |0) and 
It follows that, written as a binary expansion, qj = 0.ajv_i . . . ag = ajv_i2^^ + - ■ -+002"^. 
There is no unique way to quantize a classical map. Here we adopt the quantized baker's 
map introduced by Balazs and Voros which can be written as 

T = Fj^\l ^ Fn-i) , (5) 

where F^_i acts on the — 1 least significant qubits, and / is the identity operator acting 
on the most significant qubit. Saraceno has introduced a quantum baker's map with 
stronger symmetry properties by using antiperiodic boundary conditions, but in this article 
we restrict the discussion to periodic boundary conditions as used in . It is straightforward 
to adapt the discussion in this paper to Saraceno's version of the map P^ . 
Note that |2|] 

\^) = Fj^^\aN-i) ® ■ ■ ■ ® |ao) 

= (|0) + e-2-^(0-'*")|l)) ® . . . ® (|0) + e-2-^(0-«iV-2...ao)|l^) ^ (|Q^ + g-27r*(0.a^_l...ao)|^^) ^g) 

and 

10) = {I® FAr_i)"Vw-l) ® • • ■ ® ko) 

= \aN-i) ® (|0) + e~2-(0-'^o)|l)) ® ■ ■ ■ ® (|0) + e^^^^C^-'v-^.-ao)!!^) _ (7) 

As lip) = T\(j)), T can be seen to perform a shift of the qubits ||^, in which the most 
significant qubit of the argument |0) is transformed in a way that depends on all the other 
qubits, becoming the least significant qubit of the image The quantum baker's map 
is thus equivalent to a shift map on a quantum spin chain ^j, in analogy to the symbolic 
dynamics for the classical baker's map. 

The quantum baker's map can be realized using the following basic unitary operations 
or quantum gates: the Hadamard transform gate acting on the mth qubit and defined 
in the basis {|0), |1)} by the matrix 

and the phase gate BmniO) operating on the mth and nth qubits and defined by 

B:mn{.0) |ai_i) (g) ■ ■ ■ ® |ao) = e*'^'"""" \aL-i) ® ■ ■ ■ ® |ao) , (9) 
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where 



/ J if (3^ 1 , /'in^ 

(Pa^ar. ~ I Otherwise. ^ ^ 

In addition we define the gate Smn which swaps the qubits m and n. 

In D = 8 = 2^ dimensional Hilbert space, one iteration of the quantum baker's map is 
performed by the sequence of gates 

T = So2AoBU'^/2)BU^/A)A^Bl,{n/2)A2SoiAoBoi{n/2)A, . (11) 

The corresponding pulse sequence on the NMR computer is quite long and complicated 
(see Section |T| and the Appendix). Therefore we introduce a simplified version, Tj^, of the 
quantum baker's map Tm maps each of the states 

lajv-l) ® (|0) + e-2-(0.ao...a^„2)|^^) ^ QQ^ ^ g-2^i(0.ai...a^_2)|l^) 

(|0) +e-2-*(0-'^iv-2)|i^^ (12) 

to 

QQ^ _|_ g-27ri(0.ao...ajv_2a]V-i)|]_^^ (g, Qg^ _j_ g-27ri(0.ai...a]v_2ajv-l) | 

®---®(|0) + e-2'^'(°-"'^-i)|l)) , (13) 

and is thus equivalent to — 1 aAr_i-controlled rotations, a Hadamard transform on the 
most significant qubit, and then a cyclic shift of the qubits. In D = 8 = 2^ dimensional 
Hilbert space, one iteration of the map Tj^ is performed by the much shorter sequence of 
gates 

Tm = SoiSo2AoBU'rr/4)Bl{n/2) . (14) 

Like the quantum baker's map T, the simplified map Tj^ is a shift on a string of qubits, 
although Tjii leads to different phase relations between the qubits. The two maps can 
thus be expected to have similar chaotic behavior. We have confirmed this expectation by 
comparing the numerical results of Section |rV| with simulations of the full quantum baker's 
map T these simulations are not included here because, unlike the results of the present 
paper, they are based on unrealistic assumptions for the experimental parameters. 



III. IMPLEMENTATION ON AN NMR QUANTUM COMPUTER 

A. The system and its Hamiltonian 

We choose for our physical system the molecule trichloroethylene (Fig. 1), in which the 
nuclear spins of the hydrogen and two carbons serve as our qubits. These spins weakly 
interact with each other on a single molecule, but are effectively shielded from the environ- 
ment by rapid tumbling. The molecules are placed in a strong, uniform magnetic field and 
subjected to RF pulses at various frequencies. 
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We denote by X, Y, and Z the ax, Cy, and cr^ Pauli matrices, respectively, and indicate 
with a subscript (e.g., Xh) to which spin they apply. The Hamiltonian of the three spins in 
the interaction picture is 

H = ^-^ZhZc, + + Yc.Yc, + Zc.Zc,) + j^h^c, + (15) 



where ji ^ 203 Hz, j2 ~ 102 Hz, ^ 10 Hz and 5 ^ -905 Hz |25|. 

It is convenient to approximate the interaction term between the two carbons by 
(j2/4)2'c^Zc2, omitting the XX and YY terms. This approximation is somewhat diffi- 
cult to justify, but greatly simplifies the description of the quantum gates. Generally, it 
works well if the spin precession frequency of two spins differs by an amount large compared 



to the spin coupling j [26|. This is true between H and Ci (where the approximation has 
already been made), but only borderline between Ci and C2. The approximate Hamiltonian 
is then 

H' = ^-^ZhZc, + ^-fZc^Zc, + ^-fZuZc, + ^Zc„ (16) 

For our numerical simulations, we will assume the form H' of the Hamiltonian, but we have 
checked the dependence of our results on this approximation. Whenever the data curves 
obtained using H differ from those obtained using H' , we present both curves. 

Since js is small compared to the other terms in the Hamiltonian, it may be safely 
neglected for the design of the pulse sequences (see below). The js term cannot be neglected, 
however, in the simulation of the full dynamics including noise [see Eq. (0)], since it is of 
the same order of magnitude as the leading noise terms. Thus, for the discussion of the pulse 
sequences only, we will assume the following, further simplified form of the Hamiltonian, 

H" = ^-^ZhZc, + l^c. Zc, + ^-Zc,. (17) 

In addition to these interactions, we can apply RF pulses which rotate the nuclear spins 
about the x and y axes. By controlling the pulse frequencies, these can be selectively 
applied to single spins (soft pulses), or to two or three spins at once (hard pulses) 0. In 
our simulations, we use instantaneous pulses, i.e., we assume that the duration of the pulses 
is very short compared to the timescale of the Hamiltonian H" . This assumption is only 
marginally satisfied for soft pulses. The general conclusions of this paper, however, are 
not affected by this approximation, since they do not depend on the precise form of the 
implemented map. 

For the purposes of our simulations we also assume that the RF pulses are timed with 
perfect accuracy. Unfortunately, this is not the case in experiment, where one can expect 
errors of 1-10% or even higher ||2^. This is an additional complication, which muddies the 
argument without changing its basic conclusions; hence we neglect it. 

The general form of a quantum algorithm in an NMR computer is a sequence of pulses, 
causing rotations of the individual bits, interspersed by precisely-timed delays during which 
the undriven Hamiltonian couples the neighboring spins 0-§| . In describing such a sequence, 
X{9) denotes a rotation about the x axis by an angle 9. This is equivalent to multiplying 
the state by the operator exp(i^X/2). Y{9) denotes a similar rotation about the y axis. 
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U{t) indicates a delay of duration t, during which the Hamiltonian H" acts. A sequence is 
to be read from right to left, i.e., the rightmost operation is performed first. In this way, 
composition of the operations follows the same sense as operator multiplication. Subscripts 
indicate which spin is acted on. 

The basic gates which form the algorithm are constructed from these simple pulse se- 
quences, as we will show below. Note that in all cases we neglect the overall phase of the 
state. Thus, two gates will be considered equivalent if they agree up to an overall phase. 



B. One-bit gates 

We are already equipped with two families of one-bit gates, the x and y rotations. Two 
other useful gates, however, are lacking: z rotations and Hadamard transforms. Fortunately, 
in both cases these can be built from sequences of simple pulses | |28| ] . 



Rotations about z can be constructed from a sequence of three x and y rotations p9[] : 



Z{e) = X{-n/2)Y{e)X{Tr/2). (18) 

In fact, there are several such combinations which can be used to produce a z rotation. As 
gates are put together in a quantum algorithm, it is often useful to choose their precise form 
so that a certain number of pulses combine or even cancel out, hence simplifying the overall 
sequence. Thus, we can write any of the following: 

Z{e) = X {-it/ 2)Y (6) X {71/2) 
= X{n/2)Y{-e)X{-Ti/2) 
= Y{n/2)X{9)Y{-n/2) 

= Y{-7r/2)X{-e)Y{n/2), (19) 

as convenient in constructing the algorithm. 

We also need the Hadamard transform A defined in (^, which can be effected by a pair 
of X and y rotations 

A = Y{n/2)X{n) 

= X(-7r)r(-7r/2). (20) 

Again, the form is chosen to simplify the sequence as much as possible. 



C. Two-bit gates 

Given the large family of one-bit gates to choose from, we need only a limited selection 
of two-bit gates. This is good, because we have only a limited selection to choose from. 
To build our algorithm, all we need is the phase gate Bij{6) defined in equation (H) which 
couples pairs of bits. The phase gate between bits i and j can be decomposed into 

Bij{e) = exp{ieZi/A) exp{ieZj/4:) expiieZiZj/A). (21) 
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The ZiZj term is the critical one. Transformations of that type are produced by the Hamil- 
tonian evolution ([T7|). However, this includes unwanted additional terms. We can effectively 
eliminate these terms by the technique of refocusing, in which we "undo" the evolution of 
all but the selected term in the Hamiltonian pD|,^J5[]. Because of the anticommutation of 
the Pauli matrices, 

Xiexp{i9ZiZj) = exp{—i9ZiZj)Xi, i ^ j. (22) 

Thus, if we stick an x pulse in the middle of a period of Hamiltonian evolution, it can 
effectively remove the unwanted terms. 

Xh Gxp{—iH"T)XH Gxp{—iH"T) = exp(— irj22'ci^C2/2 — irSZci) 
Xc, exp{-iH"T)Xc, exp{-tH"T) = exp{-iTj,ZHZcj2). (23) 

We use this to build phase gates between neighboring spins. 

BcM-0) = Zu{-9I2)Zc,{~9I2)Xc,{j^)U{t)XcMV{t), (24) 

where r = 9/2ji. Note that one can use y rotations instead of x, and shuffle the order of 
these operations; also, one can choose any of the expressions (0) for the z rotations. A very 
similar expression applies for -BciC2! using X^ instead of X^^, but we must also include an 
additional Zc2 rotation to undo the effects of the bZc^l^ term in the Hamiltonian. 

While the phase gate is useful in itself, it can also be used to produce the controUed-not 
(CNOT) gate by nesting it between two Hadamard transforms, Cij = AiBij{n)Ai. 



D. Swaps 

A disadvantage of the Hamiltonian (|l^ is that it only couples neighboring spins. If we 
wish to perform a two-bit gate on H and C2, we have to swap one of them with the central 
spin, Ci, using the swap gates Sc\h and 5'ciC2- 

Swaps can be built from a sequence of CNOT gates: 

= CijCjiCij. (25) 

Each CNOT is composed of a phase gate and two Hadamard transforms, as shown above, and 
the phase gates in turn are built out of precisely timed Hamiltonian evolution interspersed 
with RF pulses. Thus, to swap the H and Ci spins requires the sequence 

= ^H^i/(7r/2)Zc,(7r/2)Xc,t/(r)Xc,t/(r)AHAc,ZH(7r/2)Zc,(7r/2)Xc,t/(r)Xc,t/(r) 
xAHAc,ZH{7r/2)ZcA^/2)Xc,U{T)Xc,U{T)AH, (26) 

where r = tc /2ji. The C1-C2 swap is similar. Since each of the Hadamard transforms and 
z rotations is itself a product of several RF pulses, we see that the swap gate is quite large 
and complicated. What is more, the three phase gates take considerable time, allowing the 
system to be affected by decoherence. The number of RF pulses can be somewhat reduced 
by carefully choosing the expressions for the Hadamard and z gates to cancel as many pulses 
as possible, but there is nothing to be done about the time for the phase gates. 

It is clearly to our advantage to perform as few swaps as possible. We follow this precept 
in designing the map T^vi- As it is, more than half the time of the algorithm is spent in 
swapping. 
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E. The map Tm 



In section II we defined a simplified version of the quantum baker's map, which in the 
three-bit case is given by (^^. This map has two important traits. First, all of its two- 
bit gates involve bit 0. Second, the two swap gates perform a cyclic shift of the bits, 
^ 2 ^ 1 ^ 0. 

Since we can only couple neighboring spins, it makes the most sense to make the central 
bit, Ci- We can then identify bit 1 with H and bit 2 with C2. In terms of the physical bits, 
then, the map becomes 

Tm = Sc,HSc,c,Ac,Bc,c,i-Tr/A)BcM-^/'^)- (27) 

Further examination reveals that the SciC2 is unnecessary. The shift 1 ^ is necessary in 
order to keep bit in the central position. The labeling of bits 1 and 2, however, is arbitrary. 
We could just have easily identified bit 1 with C2 and bit 2 with H. Thus, we need not 
perform an actual physical swap gate here; a mental re-labeling of the bits is sufficient. 

This means that the labeling of the bits will switch back and forth between even and 
odd iterations of the map. Thus, the physical realization of the map will also differ between 
even and odd steps. We therefore have two alternating sequences, 

Todd = Sc,hAc,Bc,c,{-7^/^)BcM-^/'^) 

Teven = Sc,C,Ac,Bc,H {-^^ / ^) Bc,C,i-^ f^) . (28) 

Each of these gates represents a given series of RF pulses and delays. By choosing the form 
of the z rotations and A gates carefully, and ordering the operators for the B gates appro- 
priately, a certain amount of cancellation is possible, simplifying the sequence somewhat. In 
this way we arrive at the following sequence of elementary pulses: 

Todd = XH(-37r/2)r^^(-7r/2)Fc,(7r/2)Xc,(-7r/2)Fci(-vr/2)t/(ri)Xc,(7r)t/(ri)X^,(-37r/2) 
xXcA-^n/2)YHi-n/2)YcA-^/2)Uin)Xc,in)Uin)XHi-Sn/2)XcA-^n/2) 

xYHi-7T/2)YcA-^/2)Uin)Xc,{n)U{n)Xc,{7r/2)Yc,i5n - tt/S) 
xXc,i7r/2)XH{-57r/A)YHi-n/2)Xc,i-Un/8)Yc,i-n/2)Uin), (29) 

where Ti = 7r/2ji. Note that because ji ~ 2j2 we can combine the two B gates into a single 
time delay for Todd- 

This unfortunately works in exactly the wrong direction for Tcven? which consequently 
makes the sequence somewhat longer and more complex: 

Teven = Xc,{26t2 - St, /2)YaA7r /2)Yc,{7i /2)Xc,{7c /2)Yc,{7c /2)U {r2)X h (tc) 

xUiT2)XcA-^7T/2)Xc,i26r2 - 3Tr/2)Yc,i-7T/2)Yc,i-n/2)Uir2)XcA^)U{T2) 
xXc,{-37f/2)Xc,{25r2 - 3n/2)Yc,i-n/2)Yc,i-n/2)U{T2)Xc,in)Uir2) 
xXHin/2)YHi-n/8)XHi-7r/2)Xc,{ASrs - / A)Yc, {- n / 2) Xci-lln/ 8) 
xYc,{-n/2)U{3T,/2)XH{n)U{5vs/2), (30) 

where T2 = 2ti and r^, = ti/2. (Recall that these sequences should be read from right to 
left.) The total delay time for Todd is 7ri, while the delay time for Tgven is 14ri, taking twice 
as long. Thus, intrinsic noise affects the system more during even steps than odd ones. 
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F. The regular map T-ji 



Having defined a quantized version of a classically chaotic map, we need to define a 
regular map for comparison. This map should take the same time per iteration as the map 
Ta4, on average, so that it is affected equally by the intrinsic noise. 

The simplest possibility would be a map that does essentially nothing. This implies that 
the effects of the Hamiltonian in coupling the spins must be suppressed. We do this by 
refocusing, as described above. Consider the sequence 

Tn = Xc,(7r)f/(r4)Xc,(7r)t/(r4)Xc;,(7r)f/(r4)Xc,(7r)t/(r4)Xc,(7r)f/(r4)Xc,(7r)f/(r4) x 

XcA'^)U{r,)XcMU{n). (31) 

The X rotations flip the sign of all terms containing in the Hamiltonian, so that the 
effects of these terms cancel out. That is, all but the dZc^/'^ term of ( ]T7| ) are effectively 
eliminated: 

Xc,(vr)f/(r)Xc.,(vr)f/(r) ^ exp(-z5rZcJ. (32) 

Thus, we see that 

T-ji K. exp(-4z(5r4Zc2). (33) 

By performing an appropriate z rotation at the end we could eliminate this cumulative 
rotation of the C2 spin, but as it makes no difference to the regularity of the map we leave 
it unaltered for simplicity. The total time 8x4 should equal the average iteration time of the 
map Ta4, so = 21ri/16 = 2l7r/32ji. 



G. Noise and the master equation 

While the nuclear spins are fairly well isolated, this isolation is not perfect. Interactions 
with the external environment cause the spins to deviate from strict Hamiltonian evolution, 
a deviation which becomes increasingly important as the length of the NMR calculation 
increases. 



The effects of the environment are described by two timescales Ti and T2 Ti is 

the timescale on which the spins relax to the thermal state. T2 is the timescale of phase 
decoherence, the result of the spins becoming correlated with the state of the environment 
and hence going from a superposition of states to a statistical mixture. For the molecule 
trichloroethylene, Ti exceeds T2 by an order of magnitude, and can be safely neglected for 
this problem. 

We model the effect of T2 by replacing the state vector \ip) with a density matrix p and 
the Schrodinger equation with a Markovian master equation in Lindblad form ||11|| , 



^ = -t[H\ P] + ^ E 2 {LkpLl - LlL.p - pLlL,) , (34) 

where H' is the Hamiltonian (^) and the are a set of operators chosen to model the 
effects of decoherence. The simplest choice of L^'s which capture the essential physics are 
proportional to the Z operator for each of the three spins: 
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Li = ^|T„Z 
L3 = x/^CiZci- 



(35) 



The decay rate F is proportional to I/T2 for each spin. 

Using quantum trajectory techniques PT] , |52| , we can solve the above master equation 
numerically for various choices of T2 and various initial states. We treat the RF pulses as 
instantaneous unitary transformations, just as before. They interrupt the continuous master 
equation evolution given above. 



IV. MEASURES OF CHAOS 

A. Rate of entropy increase 

The Kolmogorov-Sinai (KS) entropy, equal to 1 bit per step for the classical baker's 
map, measures the asymptotic rate at which information about an initial phase-space point 
must be supplied in order to be keep the ability to predict the nth iterate of the map to 
a given accuracy. Furthermore, if a stochastic perturbation is added to the map, the KS 
entropy measures the average entropy increase per step that results from averaging over 
the perturbation. Both of these properties are very easily understood in the shift-map 



representation of the map |20 



For quantum systems, a constant rate of entropy increase in the presence of environment- 



induced decoherence has been proposed |T2|JT3[] as a signature of chaos. This rate of entropy 



increase is closely related to the concept of quantum dynamical entropy [^,14|. As in the 



classical case, the analysis is greatly facilitated by the shift-map character of the quantum 
baker's map. 

The entropy increase in a chaotic system is due to the exponential magnification of 
small-scale, local perturbations. Since perturbing all qubits in the spin chain at the same 
rate corresponds to a perturbation on all scales simultaneously, an entropy increase can be 
regarded as a signature of chaos only if the qubits are perturbed at different rates. This 
condition is fulfilled in our system, since the spin relaxation timescales T2 for the hydrogen 
and carbon atoms in trichloroethylene differ by about one order of magnitude. 

In our first proposed experiment, the system is prepared in the initial pure state 

1^0) = |0,) ® |0,) ® |0,) = -1=(|0) + ® (|0) + ® (|0) + (36) 

To find the system entropy, S{n), after n steps (n = 1, . . . , 6), the map Tm is iterated n times 
by applying the pulse sequences ( PP| ) and (^) alternately. The final density operator p{n) 
is then measured using quantum tomography and the entropy is determined from 

S{n) = —tr [p(n) log2p(n)]. Although standard NMR techniques give one the traceless part 
of the density operator only, there exist methods (at least in principle) to determine p(?7,) 
fully, as is required for determining the entropy S{n). One such method is to do tomography 



of the quantum operation [0], where the pulse sequence is applied to different initial states. 



10 



Another possibility would be to estimate the relative size of the traceless component by 
measuring the signal strength. 

We have simulated this experiment by numerically solving the master equation (0). The 
result of the simulation is shown by the data points labeled "chaotic" in Fig. 2. The data 
points labeled "regular" are the results of a control simulation using the regular map T-]i 



defined in Section [111 F| . In both cases, the entropy increases rapidly and approaches the 
value of 3 bits, which is the maximal von Neumann entropy in 8-dimensional Hilbert space. 
There is no clear difference between the regular and chaotic cases. The apparent reason is 
that due to the small T2 timescales l/Fcj and 1/Tc2 two of the three qubits are strongly 
perturbed. To show that our conclusions do not depend on the assumptions leading to the 
simplified Hamiltonian (|16D, we have repeated the simulations leading to Fig. 2 using the 
full Hamiltonian (]T5|). The data points labeled "regular +xy" show the result in the regular 
case; in the chaotic case, the extra terms in the Hamiltonian made no significant difference. 

If only one qubit is strongly perturbed, however, there is a clear difference between 
chaotic and regular behavior. In the simulation shown in Fig. 3, both l/Fc^ and l/F// are 
relatively large compared to the total delay times needed for the sequences Todd and Teven 
(see Sec. piIE| ), i.e., only the C2 spin is strongly perturbed. These T2 timescales cannot be 
achieved with the molecule trichloroethylene, but there may exist other molecules with the 
desired properties. In this case the "regular +xy" plot differs more than in Fig. 2, but is still 
clearly distinguishable from the chaotic case. 

In Fig. 3, the entropy increase in the chaotic case does not have a well-defined linear 
regime. The reason for this is the relabeling of the qubits at each step, which was introduced 
in Sec. [HI E| to reduce the complexity of the pulse sequence. At alternate steps, the strongly 
perturbed C2 spin thus represents either qubit 1 or qubit 2. One could eliminate this effect 
by performing extra physical swap operations as described in Sec. [HI E[ Another possibility 
is to introduce artificial perturbations. 

One can apply an artificial perturbation to a map T by adding an extra z rotation to the 
least significant bit, producing the perturbed map T' = e'^^^/'^T. At each step, one randomly 
chooses either the perturbed or unperturbed map. A density operator results from averaging 
over the two possible outcomes. This is equivalent to applying a superoperator P, 

V{p) = ^(p + e'^^'/^pe-'^^''^). (37) 

The results of this and the next section depend on the strength and locality of the 
perturbation, but are rather insensitive to its exact form. This is a general property of 
entropy measures for quantum chaos, which have to be defined with respect to a class of local 
perturbations (see, e.g., P). Clearly, a perturbation that commutes with the unperturbed 
dynamics will not reveal any chaotic properties of the latter. 

To perturb the same logical qubit at even and odd steps of the map Ta^, we apply after 
each odd step the perturbation superoperator 

Podd(p) = \{p + e^-^-'^pe-^^^-/') (38) 
and after each even step the perturbation superoperator 

Peven(p) = \{p + e^^^C, /2^g-..Ze, /2) (39) 
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to the density operator p. For the regular map Tt^, we apply Vodd at each step. In an actual 
experiment, a convenient way of averaging over different perturbations consists in applying 
selected gradient fields . 



The results are shown in Fig. 4. We have assumed large relaxation times for all three 
spins. The extra terms in the Hamiltonian had only a slight effect on the chaotic case, and a 
negligible effect on the regular case. The simulation differentiates well between the regular 
and chaotic cases, and the latter shows the expected linear increase in entropy, followed by 
saturation at 3 bits. Unfortunately, this simulation assumes possibly unrealistic relaxation 
times. 



B. Hypersensitivity to perturbation 

Hypersensitivity to perturbation is an information-theoretic criterion for classical and 



quantum chaos |15|-[17| which has been shown to be equivalent to a standard definition of 
classical chaos under general assumptions |T^. Suppose a system is perturbed, for instance 
by being acted on by an unknown force with a known distribution. Averaging over all 
possible perturbations causes the entropy of the state to increase. One can reduce this 
entropy growth by obtaining information about the perturbation, such as the actual value 
of the force to some precision. By having more information about the perturbation, the 
uncertainty in the state (and hence its entropy) decreases. To reduce the entropy growth by 
an average amount AS requires information / > AS about the perturbation. In particular, 
we want to know the minimum amount of information Jmin needed to produce a given entropy 
reduction AS. A system is hypersensitive to perturbation if the information Jmin needed to 
lower the system entropy increase by an average amount AS is very large compared to AS*. 



Precise definitions of the quantities Jmin and AS* are given in [0]. For a general introduction, 
see 



We will show in this section that hypersensitivity to perturbation can be detected in the 
3-qubit quantum baker's map even in the presence of the actual noise levels for trichloroethy- 
lene. As in the last section, we define perturbed maps = e^'^^^'^T-^ and Tj^ = e^'^^^^'^T^^^ 
or T'j^ = e'^'^^'^^/'^Teven for odd or even steps, respectively. Applying randomly at each step 
either the perturbed or the unperturbed map leads, after n steps, to 2" possible different 
perturbation histories. 

Due to the fast decoherence, it is necessary to limit the number of steps to n = 3, 
corresponding to A^hist = 2" = 8 different perturbation histories. The proposed experiment 
is to apply, each time starting from the initial state iV'o) defined in (|36D, all Ahist = 8 
perturbation histories, to obtain the list of final density operators C = {pi, . . . , PNi^i^t} by 
quantum tomography, and to analyze the distribution of the Ahist density operators in 
density operator space. We assume that, in a random trial, all Ahist perturbation histories 
would occur with the same probability 1/Ahist- We can find the entropy S^^ax = — ti'(plog2 p) 
of the average density operator 



P = 1^ E Pr (40) 



A^hist ,-=1 



As argued in the last section, this should grow quickly with the number of iterations n. 
Our simulation using the relaxation times for trichloroethylene (Fig. 5) gave, for n = 3, 
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'S'max = 2.67 bits in the chaotic case, S'max = 2.74 bits in the regular case, and S'max = 2.72 
bits in the regular case with x and y terms: nearly identical values. The slightly lower 
entropy in the chaotic case arises because of the alternation between Todd and Teven! it is 
otherwise of no significance. 

The values for the entropy increase alone, therefore, do not reveal much about the dis- 
tribution of the density operators in the ensemble. In particular, they do not reveal whether 
the ensemble is orthogonal (in which case the entropy increase corresponds to purely classical 
information) or nonorthogonal (corresponding to quantum information). Obtaining infor- 
mation about which perturbation history has been realized can reduce the entropy from 
•S'max to a lower value S; analyzing the dependence of /S.S = S'max — S on the information 



needed gives a measure of how nonorthogonal the ensemble is |T7|. 

We could obtain the total possible information by determining exactly which pertur- 
bation history occurred. This corresponds to logg A^hist = 3 bits of acquired information, 
and would reduce us from considering the average density operator p to considering only a 
single final density operator pj. However, we could also obtain pariia/ information about the 
perturbation history by partitioning the A^hist final density operators into R < A^hist groups, 
and determining only which group the operator was in. Since we are actually interested in 
the minimum information Jmin needed to produce a given entropy reduction, we would like 
to choose groupings which maximize the entropy reduction. 

More precisely, consider a partitioning of the list C into R groups, labeled hj r = 1, . . . , R. 
We denote by A^^ the number of density operators in the rth group {Y.r=i = A^hist)- The 
A^r- density operators in the rth group and their probabilities are denoted by p^, . . . , p^^, and 
ql, . . . ,g]v^, respectively. In our case, all g[ = 1/A^hist- In a random trial, the system state 
will be in the rth group with probabihty 



Nr 



Pr=T.<ll = W-- (41) 

The knowledge that the system state is in group r is described by the density operator 

Nr Nr 

i=l i=l 



We define JT^ the system entropy conditional on being in group r, 

Sr = -tr{pr log2 pr) , (43) 

the average conditional entropy 

S = ^PvSr = S'max - AS < S'max , (44) 
r 

and the average information 

I=-Y.Pr^Og,Pr>AS. (45) 

r 

The information Jmin needed about the perturbation to reduce the system entropy by an 
amount AS' is now defined as the minimum of / over all groupings for which AS' > AS, i.e., 
all groupings for which the system entropy is reduced by at least AS. 
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The particular case we are treating is simple enough that we could actually try all possible 
groupings to find the one which minimizes / for a given AS". In general, however, for large 
numbers of iterations the number of groupings grows far too rapidly to exhaustively consider 
all possibilities. Instead, we must find an efficient grouping algorithm which approximates 
this minimum. 

To find an approximation to /min as a function of AS*, we introduce the concept of nearly 
optimal groupings. Given a tolerable entropy AS, we want to partition the list of density 
operators C into groups so as to minimize the information / without violating the condition 
AS > AS or S < S'max — AS. To minimize /, it is favorable to make the groups as large as 
possible. Furthermore, to reduce the contribution to 5 of a group containing a given number 
of density operators, it is favorable to choose density operators that are as close together as 
possible in some suitable sense (see below). 

To find a nearly optimal grouping into R groups, we first choose R density operators 
at random from the list C. Then for each of the remaining density operators in the list, 
we execute the following procedure. Let pk be the next density operator in the list to be 
grouped, and let denote the average of all density operators grouped into group i so far 
(i.e. excluding all those that have not yet been grouped). Then p^ is added to that group j 
for which the "distance" 

1 1 1 r ] 

dip'j, Pk) = - - tr-(p^. + Pk) log -(Pj- + Pfc) - - 2 [ - -trpfc logpfc - -trp^. logp^J (46) 

is minimal. Of course there exist many alternative grouping algorithms, of which we tried 
several, but the one described above gave consistently the best results (i.e., the smallest / 
for a given AS"). 

Figure 5 shows /min versus AS for both the chaotic and the regular case. The slope of the 
chaotic curve is roughly equal to 6, i.e., about 6n bits of information about the perturbation 
are needed to reduce the system entropy increase by n bits. In the regular case, 1 bit of 
information about the perturbation is sufficient to reduce the system entropy increase by 
almost 0.7 bits (0.5 bits with XX and YY terms). The criterion of hypersensitivity to 
perturbation thus differentiates well between chaotic and regular behavior. Furthermore, 
the slope of 6 in the chaotic case is not very far from the dimension of Hilbert space, D = 8. 
A slope close to D is characteristic for a random distribution of pure states in Hilbert space, 
and has been conjectured to hold for chaotic quantum systems [^. The steep slope indicates 
that the ensemble is highly nonorthogonal. 



V. CONCLUSION 



The quantum baker's map can be implemented with present-day technology on a 3- 
qubit NMR quantum computer. In order to investigate the feasibility of quantum chaos 
experiments using this system, we have numerically solved the master equation for the 
NMR system, including the Hamiltonian time evolution, the RF pulses, and phase noise due 
to the environment. 

We have proposed and analyzed two specific quantum chaos experiments. In both exper- 
iments, we compare the quantum baker's map with a trivial map. One experiment analyses 
the increase of the von Neumann entropy due to decoherence. We show that in principle 
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this experiment distinguishes well between the chaotic and regular cases, but a successful 
execution requires lower decoherence rates than seem to be achievable at present. 

The second proposed experiment looks for hypersensitivity to perturbation, an 
information-theoretic criterion for chaos. We have shown that hypersensitivity to perturba- 
tion can be detected in the 3-qubit quantum baker's map even in the presence of the actual 
noise levels for trichloroethylene. Using realistic estimates for the experimental parameters, 
our simulations show that this criterion differentiates very well between chaotic and regular 
behavior. 

We have thus shown that the quantum baker's map displays behavior of fundamental 
interest even for the 8-dimensional Hilbert space of three qubits. Quantum computers can 
be used to study quantum chaos under highly controlled experimental conditions. 
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APPENDIX A: 

Here we give the RF pulses corresponding to the baker's map T defined by the sequence of 
gates (0). This unfortunately includes interactions between non-neighboring bits. However, 
one can get around this problem by inserting an extra pair of swap gates 5*01 at the beginning 
and end of the iteration, making the gate sequence 

T = S'oi5i2Ai4i(7r/2)42(vr/4)AoSoi42(vr/2)A2Ai5oi(7r/2)AoS'oi. (Al) 

The 5*01 gates cancel between iterations, so one need only swap at the beginning and end of 
the entire run; and since the labeling of bits is arbitrary, these swaps can be absorbed into 
the process of initial state preparation and final state tomography. Thus, for each iteration 
we perform the sequence of gates 

r = 5i2Ai4i(7r/2)Sj2(vr/4)/lo5oi5l2(vr/2)A2AiSoi(7r/2)Ao, (A2) 

which only couples neighboring bits. 

Here there is no need to change representation every other step, so we may fix a label 
onto the three bits. In this case, we will identify bit with if, bit 1 with Ci, and bit 2 with 
C2. The sequence of pulses corresponding to each gate is described in detail in section HI, 
above. All that is required is to combine them into a pulse sequence for the entire map: 

r = YcA'^/2)XH{Tr)Xc,{Tx)U{AT)XH{iT)U{AT)YcA^l2)XcA^r6)YcA^r6 - tx/2) 
xXc,(7r/2)Xc,(7r/2)f/(4r)X^,(7r)?7(4r)rci(7r/2)rc.(^/2)Xc;,(7r)Ac,(7r)[/(4r) 
xXh {tx)U{At)Xc, (-7r/2)Xc, (-7r/2)Fci (-3vr/4)Xc;, (7r/2)Fc. (85r - ^/2)Xc, (-7r/2) 
xf/(r)Ac,(7r)f/(r)Ac,(7r/2)A^(7r/2)r^(7r/4)A^(7r/2)Fc,(7r/8)Xc,(-7r/2)f/(r) 
xA^(7r)f/(r)Xc;,(7r/2)rc;.(vr/8 - 25r)Xc;,(7r/2)[/(2r)Xc,(7r)[/(2r)Xc;,(7r/2) 
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xXj^(7r/2)C/(2r)Xc,(7r)C/(2r)XH(-37r/2)y^^(-7r/2)Xc,(-37r/2)yci(-7r/2)C/(2T) 
xXc,{'K)U{2T)YH{Tr/2)U{2T)XH{'K)U{2T)YcAT^/2)XcM5T - 7r/4)Xc,(-7r/2) 
xyci(-7r/4)Xci(-57r/4)yci(-7r/2)C/(3T)Xc,(7r)C/(3T)y^(7r/2)X^^(7r/4), (A3) 

where the basic timescale is r = 7r/4ji = ri/2. 
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FIGURES 



FIG. 1. The molecule tricliloroethylene. 

FIG. 2. Entropy versus number of steps for the regular map T-pi and the simplified baker's 
map Tm- The decoherence times are l/Fc^ = 0.7s, 1/Th = 4.0s, and l/Fc^ = 0.4s, i.e., realistic 

values. The curve labeled "regular+xy" was generated with the XX and YY terms included in 
the Hamiltonian. For the chaotic map, the efi^ect of the extra terms is negligible. 

FIG. 3. Entropy versus number of steps for the regular map T-ji and the simplified baker's 
map T_\4. The decoherence times are l/Fcj = 10s, 1/Th = 10s, and l/Fcj = 0.2s, i.e., idealized 
values. The curve labeled "regular+xy" was generated with the XX and YY terms included in 
the Hamiltonian. For the chaotic map, the effect of the extra terms is negligible. 

FIG. 4. Entropy versus number of steps for the regular map T-jz and the simplified baker's map 
Ta/( in the presence of an artificial perturbation as described in the text. The decoherence times 
are l/Tc^ = 1/Th = 1/Tc2 — 10s, i.e., idealized values. The curve labeled "ehaotie+xy" was 
generated with the XX and YY terms included in the Hamiltonian. For the regular map, the 
effect of the extra terms is negligible. 

FIG. 5. Information needed about the perturbation versus entropy reduction for the regular 

map T-ji and the simplified baker's map Tyv^. The decoherence times are l/Fcj = 0.7s, 1/Th = 4.0s, 
and 1/Tc2 = 0.4s, i.e., realistic values. The curve labeled "regular+xy" was generated with the 
XX and YY terms included in the Hamiltonian. For the chaotic map, the effect of the extra terms 
is negligible. 
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Figure 1 (Brun) 



Figure 2 (Brun) 
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Figure 4 (Brun) 
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Figure 5 (Brun) 
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